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Abstract 

We discuss the order a1 corrections to the longitudinal spin structure 
function gi(x,Q'^,m'^) which are due to heavy flavour contributions. Here 
Q denotes the virtuality of the photon and m stands for the heavy flavour 
mass. Since the exact heavy quark coefficient functions are not known yet 
we have used the asymptotic forms which are strictly speaking only valid 
in the region ^ m?. However an analysis of the exact and asymp- 
totic expressions for i^2^^*-*(x, Q^, m^) and {x , Q"^ , rn^) reveals that the 
asymptotic forms can be also used at smaller Q^. It appears that for the 
region 0.01 < x < 0.1 the NLO charm quark component can become twice 
as large as in LO. However it is still much smaller than gi{x,Q'^) due to 
light parton contributions. This is in contrast to the observations made for 
F2{x,Q'^) at small x. Also the charm quark contribution to the Bj0rken 
sum rule turns out to be very small. 



1 Introduction 

Heavy quark production in deep inelastic electron-proton scattering proceeds via 
the following reaction (see Fig. |^) 

e-{h) + P{p) ^ e-{h) + Q(pi) (Q(pi)) +' X' . (1) 

Where 'X' denotes any inclusive hadronic final state and V stands for the neutral 
intermediate vector bosons given by 7, Z. In this paper we treat the case that 
Q'^ <^ M| so that the process in Fig. |l| is dominated by the one photon exchange 

^Talk presented at the Workshop on "Deep Inelastic Scattering off Polarized Targets: Theory 
meets Experiment" ,DESY Zeuthen, Germany, September 1-5, 1997 
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Figure 1: Kinematics of charm production in deep inelastic electron-proton scat- 
tering 



mechanism only. If we also integrate over the momentum pi of the heavy (anti) 
quark Q the unpolarised cross section is given by 



dxdy {Q^y 



■S. 



eP 



{1 + (1 - yfjF^ix, Q\ rr?) - y^Fi^{x, Q\ 



(2) 



Here S^p denotes the centre of mass energy of the electron-proton system and 
the heavy quark component of the spin averaged structure functions are given 



by Fi{x, 



2,L). In the case the proton is polarised parallel (— >) or 



anti-parallel (■*— ) with respect to the spin of the incoming electron we obtain the 
cross section 



dx dy dx dy 



{2-y}gi{x,Q\m^) 



(3) 



where gi{x,Q'^,m'^) denotes the heavy quark component of the longitudinal spin 
structure function. Further we have defined the scaling variables 



X 



2pq 



y 



pq_ 
ph 



-Q' < 0. 



(4) 



In QCD the twist two contributions to the charm component of the structure 
functions defined above can be written as 



g2 



i,q\ 1 o 1 9 



+e 
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-z'' 



7-NS/ 
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(5) 



with 



Zmax = : — ^ which follows from s = > Am? . (6) 

Further Uf denotes the number of light flavours and ^ stands for the mass factor- 
ization scale which is put to be equal to the renormalization scale. The charges 
of the light quark k and the heavy quark Q are given by Ck and eg respectively. 
Notice that n/ is determined by all quarks which are lighter than the heavy quark 
Q and Eq. has to be understood in the fixed flavour number scheme. In Eq. 
(ID the gluon density is denoted by and the singlet (S) and non singlet 

(NS) combination of quark densities are defined by 

/'(^,/^') = T.[h{z,^^') + h{z,^^')], (7) 

k=l 



The heavy quark coefficient functions Li^k and Hi^k {i = 1,2,17; k = q,g) are 
distinguished according to the production mechanisms from which they originate. 
The quantity Lj ^ is given by the processes where the virtual photon couples to the 
light quark whereas k originates from the reactions where the virtual photon 
is attached to the heavy quark. Hence these coefficient functions in Eq. (|]) are 
multiphed by e| and eg respectively. 

Up to order the heavy quark coefficient functions are given by the follow- 
ing parton subprocesses. In order we have the virtual photon gluon fusion 
mechanism given by 

^*+g^Q + Q, (9) 

S (1) 

from which one obtains the heavy quark coefficient function H^g . Here the 
superscript denotes the order in the perturbation series. In next-to- leading order 
(NLO) one has the following processes. First one has to include the virtual gluon 
corrections to reaction (P) which have to be added to the gluon bremsstrahlung 
process 

-f* + 9^Q + Q + 9- (10) 
This leads to the coefficient function Hfg'^\ Finally we also have the subprocess 

i* + q{q)^Q + Q + q{q)- (n) 

The above reaction has two different production mechanisms. The first one is 
the Bethe-Heitler process where the virtual photon is coupled to the heavy quark 
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Figure 2: The function xH^^g'^^'^'^\x, 1) {rUc = 1.5 GeV/c) for ^ iq (lower 
solid line), Q'^ = 100 (lower dotted line), = 10^ (upper solid line), ^ 
(upper dotted line). All units are in (GeV/c)^ 
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which leads to the coefficient function H^^'^'^\ Here PS means purely singlet since 
this function originates from the partonic process where a gluon (flavour singlet) 
is exchanged in the t-channel between the light quark q and the heavy quark Q. 
The second one is the Compton process where the virtual photon is coupled to the 
light quark from which one obtains -Z^^^g'*'^'* = Lh',!^\ The expressions for the order 
as contribution H^^^ can be found in [0] {i = L, 2; unpolarised) and (z = 1; 

polarised). The exact analytic expressions for I/^g''-^^ are calculated in (i = 
L, 2; unpolarised) and (i = 1; polarised). The exact heavy quark coefficient 
functions H^q'^"^^ and Hfj^"^^ have been only calculated in the unpolarised case 
{i = L,2) (see 0) but the exact expressions in polarised scattering {i = 1) are 
not known yet. However in and |Q we could derive analytic formulae for the 
asymptotic heavy quark coefficient functions up to order for unpolarised and 
polarised electroproduction respectively. The asymptotic heavy quark coefficient 
functions are defined by 



'■i,k V"^' ^9 ' 9 ' o 



rrcxact/ 

^^i,k V"^' ^9 5 9 



(12) 



These asymptotic expressions are very useful because 

1 They provide us with a check on the exact calculation. 

2 The exact calculation of Hf^'^'^^ and Hfg"^^ {i = 2, L) is semi-analytic and 
they are only available in the form of long computer programs (see 0). 
The latter show numerical instabilties at ^ which are remedied by 
using in this region H^^™"^ . 

3 Further if these asymptotic heavy quark coefficient functions are substituted 
in Eq. (|^) it turns out that (see [§[, [0]) 

F,X'^(x,Q^m^) 

tends to unity as > 10 and x < 0.01. In the case of charm production 
(Q = c and = 1.5 GeV/c) this happens for > 20 (GeV/c)^. This 
implies that above a rather low value of it does not make any difference 
whether we use the asymptotic or the exact expressions for Hi k- 

4 The asymptotic forms of the heavy quark coefficient functions also play a 
major role in the derivation of the variable number flavour scheme (VFNS) 
from the fixed flavour number sheme (FFNS) (see 0, |@]) 

The asymptotic heavy quark coefficient functions get the typical form 
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with an analogous expression for L^^^™''. There are two ways to determine Eq. 

1 Evaluate the phase space integrals and the Feynman integrals in the usual 
way for ^ m^. 

2 Use operator product expansion techniques and the theorem of mass fac- 
torization. 

The last method was applied in and 0] to compute the asymptotic forms 
of the heavy quark coefficient functions for unpolarised and polarised scattering 
respectively. 



2 Order contributions to gi(x,Q'^^m' 

In this section we will analyse the charm component of the longitudinal spin 
structure function gi. To compute the latter in NLO the following coefficient 
functions are available: ff^^^'^*'^^) (Born), L'j^g"^''*-^-* (Compton process), H^g^^'^'^^ 

(gluon bremsstrahlung) , fj^^^^'^'^^ (Bethe-Heitler). Like in unpolarised scattering 

the gluonic coefficient functions indicated by Hi^g are the most important ones. 
Therefore deep inelastic electroproduction of charm quarks provide us with an 
excellent way to determine the gluon density in polarised scattering (see Eq. (|)). 
However as we will show below the charm component in the spin case is much 
smaller with respect to the light parton contribution to the structure function as 
observed for F2 in unpolarised scattering. This is revealed by H^g^"^'^"^^ in Fig. 
1^ where we observe a strongly oscillatory behaviour which is in contrast to its 
unpolarised analogue shown in Fig. 8 of 0. Further it satisfies the sum rule 

dxHPg{x,^,^) = 0, (15) 

which holds in any scheme. Notice that the Born contribution Hi^g also satisfies 
the above relation and we may safely assume that Eq. ( plSj) holds in all orders 
for general and m^. Further it turns out that the coefficient functions iff]. 
{k = q, g) are much smaller than their unpolarised analogues. In particular the 
large logarithmic terms of the type which are characteristic of the latter do 
not show up in the spin case. The absence of these corrections and the oscillatory 
behaviour leads to much smaller charm quark contributions to gi than observed 
for F2 or Fl. To give a better estimate of the order corrections we have to 
improve the asymptotic form of the coefficient functions ff^J^^'^"^^ and ff^^^^P'*^^^ 
a little bit. In § we multiplied them by ^(1 - ^) and added to iffJ'^P'^^) the 
soft gluon bremsstrahlung term which is universal (see Eqs. (5.1)-(5.3) in |Q). 
These approximations were tested for ^2'^^*-' and g\''^ and we found a fairly good 
agreement with the exact results as long as x < 0.1 and > 10 (GeV/c)^. 
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Figure 3: The charm component of xgi{x, Q^, ml) at = 50 (GeV/c)^. Dotted 
hne: x(7™'^*(Born); sohd hne: xgi^^^°^{NLO) 

Choosing the parton density set in [1^ (Standard Scenario) with Uf = 3 
and A = 200 MeV we have plotted the charm component of (7^^'"'^* (Born) and 
gi^^^°^(NLO) at = 50 (GeV/c)^ in Fig. 3. From this figure we infer that the 
charm component becomes maximal in the region 0.01 < x < 0.1. Further we 
observe that in NLO the result becomes twice as large than in LO. We also made 
a comparison in Fig. 4 with the quantity 5f}^^'^*(NL0) which is due to light partons 
(m, d, s and g) only. Here we see that the charm component amounts to about 
4% of (^I's^XNLO). This is in contrast to the observation made for the structure 
function F2 in unpolarised scattering (see |TT|, []T2|). In the latter case the charm 
component becomes large at very small x (i.e. x = 10~^) and it amounts to about 
25% with respect to the total structure function. In we also studied the spin 
structure functions plotted at = 10 (GeV/c)^ and = 100 (GeV/c)^. These 
plots show the same features as the ones made for = 50 (GeV/c)^ so that our 
conclusions are unaltered. 
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Figure 4: The charm component and the hght parton contribution to 
xgi{x,Q^,ml) at = 50 (GeV/c)^. Dotted hue: x^^pp'"" (charm, NLO); sohd 
hue: xgf''^\NLO) 



3 Order contributions to the Bj0rken sum 
rule 



Since we know the exact result for 



NS,(2) 



we can compute the order cor- 



rection to the Bj0rken sum rule due to heavy flavour contributions. The latter is 
given by 



dx 



grix,Q^m')-g^ix,Q'',m') 



6 ' ' y^'' 



(16) 



where the superscript (1) denotes the first moment (Mellin transform). In order 
«g we have the following properties 



^2 ^ 2 /rNS,(2)\(l) 

which shows the decoupling of the heavy quark when <^ m?. 



(17) 



rNS,(2)\(l) 



An 



41n( 



+ 



The light parton contribution to the Bj0rken sum rule is given by 



dx 



grix,Q-')-g^ix,Q'') 



1 I Ga 
6 ' Gv 



NS 



(1) 



(19) 



where Cf^ denotes the massless quark coefficient function. The first moment of 
the latter quantity has been calculated up to order in |jl3| and the order term 

in, 



has been estimated too (see e.g. 



It turns out that the charm contribution to 



the Bj0rken sum rule is even smaller than the the estimated order a], correction 



(see |]T5|). The contributions coming from the bottom and top quark are even 
smaller which can be attributed to the decoupling of heavy quarks shown in Eq. 
[T7| . Finally we have the relation 



NS, 



/i2 



+ hm 



NS, 



m 



(20) 



so that for ^ the large logarithm In ^ in the heavy quark coefficient 



function Lf^ (Eq. (|19])) can be absorbed in the light quark coefficient function 
Cfg. Hence the number of light flavours in the running coupling constant and in 
cfq on the righthand side of Eq. ( pO]) will be enhanced by one unit. 
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